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Abstract 

The ideal uniform two-dimensional (2D) Fermi and Bose gases are considered both in the ther- 
modynamic limit and the finite case. We derive May's Theorem, viz. the correspondence between 
the internal energies of the Fermi and Bose gases in the thermodynamic limit. This results in both 
gases having the same heat capacity. However, as we shall show, the thermodynamic limit is never 
truly reached in two dimensions and so it is essential to consider finite-size effects. We show in 
an elementary manner that for the finite 2D Bose gas, a pseudo-Bose-Einstein condensate forms 
at low temperatures, incompatible with May's Theorem. The two gases now have different heat 
capacities, dependent on the system size and tending to the same expression in the thermodynamic 
limit. 

PACS numbers: 03.75Fi,05.30.Fk,05.30.Jp 



I. INTRODUCTION 



Since their inception into the theoretical framework of physics, much work has been car- 
ried out on both the Fermi-Dirac as well as the Bose-Einstein distributions^. Much of 
this focuses on the behaviour of fermions and bosons at low temperatures, where the devi- 
ations in their behaviour arising from the distributions are significant and quantum effects 
dominate. The advent of comprehensive experimental techniques for lowering the temper- 
ature of systems to previously theoretical values over the last few decades has contributed 
significantly to this and revealed a wealth of unusual behaviour at extremely low tempera- 
tures, in particular the unusual phenomenon of Bose-Einstein condensation - a macroscopic 
occupation of the ground state of a gas of bosons at finite temperature. 

For most of this time, however, there has been only a cursory and passing interest in low- 
dimensional Fermi and Bose systems. This has been mainly because they have been viewed 
as little more than abstract theoretical and mathematical constructs. Moreover, it was felt 
that there was little behaviour of interest to be found in these lower dimensions. It is only 
over the last couple of decades that there has been a growing interest in the properties of 
gases obeying the Fermi and Bose distributions in lower dimensions. This has been driven 
by firstly, the growing practical importance of low- dimensional systems, such as for example 
in computers, as well as new experimental techniques (such as, for example, 2D adsorbed 
systems, in particular, 2D adsorbed helium film a 3 ' 4 ' 5 ' 6 ) that raise the possibility of being able 
to trap and study low-dimensional quantum systems. Typical systems being studied include 
superfluid and superconducting films, quantum Hall and related two-dimensional electron 
gases and low-dimensional trapped Bose gases. Further, the renewed theoretical interest in 
low-dimensional physics over the last couple of decades has led to a growing realisation that 
these systems possess equally fascinating and unique quantum behaviour of their own (see, 
for example, Ref. Q). In particular, recent observations of exotic quantum behavior, such 
as Bose-Einstein condensation (BEC) ^ 1 ^ and Fermi degeneracy^, in three-dimensional 
systems have led to a growing interest in lower- dimensional systems. 

One special theoretical interest in the 2D ideal gas is due to the fact that we can obtain 
completely analytic expressions for the chemical potential fi and other thermodynamic prop- 
erties. This arises, in part, from the energy- independence of the 2D density of states within 
the semi-classical approximation. Furthermore, in two dimensions, there is an unusual cor- 
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respondence between the Bose and Fermi cases, first noted by May^ and subsequently, 
explored further in later paper o 13 i 14 i 15 i 16 i 17 . 

In his seminal paper—, May noted the unexpected result that the thermal capacities 
of the ideal 2D Bose and Fermi gases are identical in the thermodynamic limit, i.e. their 
internal energies differ only by a constant temperature- independent amount. This only 
holds in the thermodynamic limit, or equivalently, the semi-classical approximation that 
the energy level spacing is negligible compared with the temperature T, allowing us to 
replace the discrete (quantum) structure by a continuum. Real-world systems, however, are 
always finite, typically of the order of 10 3 — 10 6 particles for trapped gases, and more for 
adsorbed systems. Moreover, as we shall see later, the chemical potential [i has a logarithmic 
dependence on the total number of particles N in the system. Thus, the thermodynamic limit 
is never truly reached in two dimensions and is an abstract theoretical construct, making it 
essential to consider finite-size effects. These may be experimentally significant, and lead to 
deviations from May's Theorem. 

In the following sections, we shall consider the behaviour of the uniform 2D ideal Bose and 
Fermi gases. In the next section, we shall briefly outline the gases in the thermodynamic 
limit and derive May's Theorem before examining the Bose and Fermi gases for a finite 
number of particles in Section IIII1 Though there is no BEC in the thermodynamic limit, 
we will demonstrate that there is a pseudo-condensation at very low temperatures for the 
finite 2D Bose gas. This is then used to show in Section HVl that the finite-size effects lead to 
differing heat capacities for the Bose and Fermi gases, and a breakdown in May's Theorem 
for gases with a finite number of particles. Further, the inadequacy of heat capacities as a 
means of determining phase transitions is discussed. 

II. THE THERMODYNAMIC LIMIT 

As the 2D ideal gas and May's Theorem have already been examined at length before in 
the thermodynamic limi t 12 i 13 i 14 i 15 i 16 i 17 , we shall simply draw on earlier approaches to provide 
a brief overview. Consider N free and identical particles of mass m in a 2D box of 'volume'^, 
with rigid-box boundary conditions. The energy values for an arbitrary particle then are 
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where k is the momentum characterising the energy states i,n x ,n y = 1,2,3... and 
{rix^Tiy^ are points in the positive quadrant of a unit square lattice. In the thermodynamic 
limit, the ground state energy tends to 0. For a finite system, however, due to the Heisenberg 
uncertainty principle, this is always greater than 0, as (n x ,n y ) = (1,1). So, for simplicity, 
we shift the energy scale, taking the ground state energy as 0, 

fc2p 2fc2 

Assuming that the system is in thermal equilibrium at temperature T and chemical potential 
/i, then the mean number of particles n (ei) at a given energy is given by 

n &) = e (e t -Jk B T ±1 - ( 3 ) 

Here and later, the upper (lower) signs are for fermions (bosons). The total number of 
particles iV and the internal energy E of the system are now given by 

oo 

# = (4) 
i=i 

oo 

E = ^ein(ei). (5) 

i=l 

where the sum is taken over the single particle energy states. In the semiclassical approxima- 
tion, we assume that the spacing between the energy levels of the particle states is negligible. 
This allows us to replace the sum with an integral over energy states using the density of 
states g(e) and the spin degeneracy a, 

POD 

f = a f{e)g(e)n (e) de, (6) 
Jo 

where / (e) is a general single-particle function of energy. However, this approximation is 
only exact in the thermodynamic limit when the system size is infinite, and at very high 
temperatures T. As we shall see in Section ITTT1 for finite systems at low T, these assumptions 
are invalid and the discreteness of the energy levels cannot be ignored. In other words, we 
will be at a temperature T below which the semiclassical approximation breaks down and 
modifications need to be made. However, it may be possible to compensate by modifying 
the semi-classical approximation to allow for finite-size effects. 

We now evaluate the various thermodynamic quantities for arbitrary spin, i.e. arbitrary 
a, for both the 2D Fermi and Bose gases in the thermodynamic limit. In 2D, g(e) and 



thereby, the Fermi energy Ep are given by, 

mA 2ixh 2 N , , 

9{6) =2^' £F = ~^A- (7) 

For comparison, we adopt Ep as a characteristic energy eq for the Bose gas also, allowing 
us to express later results in a unified form. In 2D, as g(e) is a constant, we have a simple 
integral for N, 

amA r°° de 
= 2^ J e^)I^BT ± \ ■ ^ 

This may be expressed using the Fermi-Dirac and Bose-Einstein functions^-, f v {z) and g v {z) 

respectively, 

1 r°° r v ~ x dr 

*W = rMjfi^i- (10) 

where T (z/) is the Euler gamma function. These can also be expressed using the polyloga- 
rithmic functions \X v {z f&^ , 

/„ («) = -Li, (-z) , (11) 
g v (z) = Li u (z) . (12) 

For N above, v = 1 and z is the fugacity (2 = e M//fcsT ). Solving (JHj) in terms of Eq and 
rearranging, we can now express the fugacity z and thereby, the chemical potential /1 as a 
completely analytic expression 

z FiJJ = e ^ sT = =F (l — e ±ec ' kBT ) , (13) 

where zp (zb) is the fugacity for the ideal Fermi (Bose) gas. As T — > 0, zb smoothly 
approaches 1. Denoting the fugacity zb by zb = e~ a , as zb — > 1, 

0! ( e ~ Q ) = -ln(l -e~ Q ) — * -In a. (14) 

iV diverges logarithmically as a — > 0. Thus, it has no upper bound, and there is no tempera- 
ture below which the ground state can be said to be macroscopically occupied in comparison 
to the excited states. Thus, there is no BEC in the thermodynamic limit. This is a qualifica- 
tion of the oft-paraded statement that there is no BEC for D < 2, where D is the number of 
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dimensions. However, as we shall see in Section ITTTl the logarithmic divergence implies that 
at sufficiently low T, bosons still crowd into low-lying states to give a pseudo-condensate. 
The internal energy may be similarly solved, 

E F = Nk B T f 2 ( e *c/ kB T _ y ^ (15) 

E B = Nk B T (^) g 2 (1 - e-cc/^) . (16) 



Now, using the following property of the dilogarithmic functio n 12 ! 19 Li 2 



we can easily show 



h[\-l\=g 2 {l-x)+ l -{\nx)\ (17) 



E F = E B + -Ne c . (18) 



i.e. for all T, Eb differs from Ep by a constant temperature independent energy, which does 
not contribute to the heat capacity. This gives May's Theorem^ 2 -, which states that the heat 
capacities of the ideal 2D Fermi and Bose gases are identical at all T and N, 

_dE B _dE F 

^V,Bose — gj, \ v — Qrp \v — °V, Fermi- U»J 

The above result indicates that the 2D ideal Bose gas does not undergo a phase transition, 
as any experimental measurement of the heat capacity would be the same as for a 2D Fermi 
gas, and would lack the characteristic bump that indicates the presence of BEC in the 3D 
Bose gas. However, as we shall see, this is because we have failed to consider the finite 
nature of any real systems. 



III. THE FINITE GAS 



The case of the finite ideal Bose gas in two dimensions has already been considered to 
varying degrees in earlier paper a 20 i 21 i 22 i 23 i 24 i 25 i 26 i 27 i 28 , and more recently, by Pathria 2 ^. Most 
have commented that it may be possible to have a pseudo-condensation in the Bose gas 
in two dimensions, with the transition temperature typically going as Tq ~ 1/ In N. Here, 
we present a simpler and more comprehensive analysis of the problem that not only offers 
far more insight into how the finite size of the system leads to the breakdown of May's 
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Theorem and a pseudo-condensation in the Bose gas, but further, in Section IIV1 discusses 
how the heat capacities for the Fermi and Bose gases differ for finite systems, and how such 
a pseudo-condensation may be experimentally measured. 

At T = 0, the 2D Bose gas is in its ground state, and we have the trivial case of ground 
state macroscopic occupation. Extrapolating, however, a range of finite temperatures must 
exist near T = 0, which may depend on N, where the number of particles in the ground 
state is also of O (N) - a macroscopic occupation. At very low T, this can be a significant 
fraction of N, and the semi-classical approach is now flawed as N is no longer continuous. 
Taking macroscopic occupation of the ground state as the defining characteristic of BEC, 
we consider the ground state separately analogous to standard treatments of the 3D Bose 
gas^2. For finite systems, the discrete level structure requires us to introduce a cutoff energy 
at the first excited level eq. The integral is now bounded and we can remove the logarithmic 
divergence. Thus, we can have a pseudo-BEC in a finite two-dimensional system. This must 
be associated with a breakdown of May's theorem at sufficiently low T for finite systems. 
We thus write (JBJ) for the Bose gas as 

olvtlA f°° de 

It should be noted that the above approximation is quite a simple and crude one, as it 
is only the leading term of the Euler-Maclaurin summation formula. However, it suits the 
purposes of our paper as our intention here is merely to highlight how even the simplest 
approximations show that there is a pseudo-condensation in the finite Bose gas. A more 
detailed approximation will only lead to minor modifications to the core results presented 
in the rest of this paper. Further, the leading correction term to the integral above in the 
Euler-Maclaurin expansion is of the order of i\q = n (eq) = (e( £l ~^/ fcsT ) . Taking the ratio 
Nq/Ni of the number of particles in the ground state to that in the first excited state, we 
can see 

^ « e J~Ji (21) 
Aq -/i v ' 

As the chemical potential fi for the 2D Bose gas is always negative and less than 0, the 
above ration is always greater than 1. Further, as T — > 0, \i — > and the above expression 
clearly goes to oo. Further, a more detailed evaluation of (J2T)) gives us 



This clearly rapidly goes to infinity as both the temperature T as well as the number of 
particles N increase. Moreover, for the regime in question, we are interested in very low 
temperatures where from the above equation (J22j) . the ratio N Q /Ni is clearly significantly 
large. 

Thus, we are justified in making our crude approximation and denoting the number of 
particles N by (|20|) above. Solving as before, we now have 

N _ _ amAksT ^ / _ 

and £\ prevents the logarithmic divergence. This is because the excited states cannot accom- 
modate all the particles at low T, and the excess are forced into the ground state, leading to 
a pseudo-BEC. Equivalently, changing variables e — > e + E\ and resetting the integral from 
zero, we may also express ()20|) in terms of the Bose-Einstein functions g u (z) as 

N-N = ^§^1 (z B e*™/*"»» T ) . (24) 
Thus, we now have an effective fugacity 

ZB _ Z > B = ZBe -3n^/2mAk B T 
_ e -3aTre c /iNk B T fa _ e -e c /k B T\ ^5) 

where Zb is as before. This summarizes the finite-size effects on the Bose gas. The extra 
term is dependent on the system size and tends to 1 as N — > oo. For finite iV and low T, it 
is less than 1, and so, the fugacity is modified. Thus, there is a pseudo-BEC here but none 
in the thermodynamic limit. This macroscopic occupation is a purely boson phenomenon, 
and so, the behaviour will be markedly different from that of the two-dimensional Fermi 
gas. 

For the Fermi gas, Pauli's principle forbids macroscopic occupation of any state. Thus, 
the first term is not divergent, zp (Eq. (JHJ)) still holds, and finite-size effects only appear for 
the Bose gas, where the excited states take as many particles as possible at low T, with the 
excess forced into the ground state. This occurs at the maximal value of the fugacity Zb — 1. 
Thus, we can define a characteristic temperature Tq that marks the onset of condensation 
where iVo ~ and z ~ 1, 

amAksTc 



N 



2ixh 2 



In (l — e - 3 K 2 h 2 /2mAk B To 



«bfM| (26) 



Ec V 37ra£(7 



for very low T. We may solve (ffijj) to give 



Tc k B W (4iV/37ra) ' (27) 

where, W (x) is the Lambert function— and is the principal solution for w in the equation x = 
we w . This is plotted in Fig. 1. It may be viewed as a generalisation to a logarithm. Thus, our 
result (J2*7)l is a refinement of the earlier noted result that the transition temperature behaves 
like Tc ~ 1/ In N. As iV — ► oo, the reduced characteristic temperature ksTc/ec — > 0. Thus, 
in the thermodynamic limit, we recover the conventional result that there is no BEC in two 
dimensions. The condensate fraction is 



No „ n ksT (ANk B T 
— = 1 In 



T 



(28) 



This only holds for very low temperatures T < Tc- As can be seen in Fig. 2, the condensate 
occurs at lower and lower T as the system size iV increases. In the thermodynamic limit, 
there is no BEC at all. 

Having shown the pseudo- condensation, we now move on to consider how the Fermi-Bose 
correspondence predicted by May is affected when the gases are assumed finite. 



IV. THE FINITE GAS: INTERNAL ENERGY AND HEAT CAPACITY 

The internal energies may be similarly found. The first term in (JHJ) for the Bose gas is 
divergent as it is the internal energy E gs associated with the anomalously large ground state. 
Separating this out, 

oo 

E * = E °° + E e{ e t -jLT_ r (29) 
i=l 

As we have set E gs as the zero-point energy, the first term is 0. Changing variables e — > £+£i, 
and converting to an integral, 

N f°° £ + £ l 



Eb ecJo e( e +*i-")/ fc s T -l de 

= — {k B T f g 2 (z' B ) + —NksTgi (z' B ) . (30) 

£c £c 

where z' B is defined by (l23j) . We again have an effective fugacity z' B and additional terms, 
which are negligible in the thermodynamic limit. By contrast, for the finite Fermi gas, there 



is no divergence in the ground state term, and we have the same expression for the fugacity 
(|T3*j) as before. The internal energies of the finite two-dimensional Fermi and Bose gases are 
then given by 

'k B r 



Nk B T 



ec 



/ 2 (' 



,ec/k B T 



1) 



E 



B 



Nk B T 2 \ 



3na 



k B T In (1 - z' B ) 



(31) 
(32) 



e c J' ' - 4 

Eq. (|18j) is no longer true, as we cannot simply express Ep in terms of E B as before. E B 
now has additional T-dependent terms that lead to a differing heat capacity. 

May's Theorem is thus a special limiting case for the 2D Bose and Fermi internal energies 
in the thermodynamic limit. Here, the additional factor e - 37rQ£ c/4A r fe s r tends ^ Q unity in the 
thermodynamic limit, recovering May's Theorem. The same occurs in the high T limit for a 
finite gas as above Tq, there is no macroscopic ground state occupation. Thus, above Tc, we 
may approximate the internal energies by (frSj) and (fT6*|) with differences becoming marked 
only at very low T and/or very small system size N. 

We now calculate the Bose and Fermi heat capacities by using the following polyloga- 
rithmic identity^^, 



d 

-r-Ll, 



+1 * 



dz z 
The heat capacity for the finite gas is then given by 



-Li n (z) 



(33) 



a 



V, Fermi 



Nkl 



CvBose 



2Tg 2 (1 



-e c /k B T) 



'C 



Nkl 



2Tg 2 (z' B ) + 



e c ln(l 



k B e £ d k i 



37ra/c£ 
16JV 



E C (37ra^ - ^ e -^clk B T-?,^ec^Nk B T^ 



k R T(l-z' 



8iVlnfl-^o 



(34) 



(35) 



"B- 1 V x — ^B) 

The Bose heat capacity Cy^ Bose now contains several terms dependent on both iV and T. It 
should also be noted that the last two terms in the heat capacity are very much smaller than 
the first two terms and become increasingly negligible as the system size increases. At high 
T and/or very large N, the additional terms tend to zero, i.e. as iV — > oo and/or T — > oo, 
7rctEc/Nk B T —>■ 0, and (J33j) reduces to the standard expression in the thermodynamic limit, 



V,Bose 



X e sc/k B T _ l 



a 



V.Fermi- 



(36) 
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This is one of the most important results of our paper, and shows clearly that May's theorem 
arises as the thermodynamic limit of (|34|) and (|35|). when the heat capacities of the Bose 
and Fermi gases become equal. Nevertheless, for most real-world 2D systems, above the 
characteristic temperature T@, which marks the onset of Bose-Einstein condensation in the 
two-dimensional ideal Bose gas, iraec / NksT is very small and so, even for finite gases, 
the thermal capacities of the Bose and Fermi systems are approximately equal above the 
condensation temperature. 

In Fig. 3, we plot the reduced heat capacities (Cy / ' NksT*) for the Fermi and Bose gases 
against the reduced temperature T* (where T* = ksT/ec) for varying system size N, using 
the full expression for both Cy^ose and Cy t p erm i. As the heat capacity is approximately 
linear close to T = in the thermodynamic limit, we would expect the plots for the 2D 
Fermi gas to tend to a horizontal line as T — > 0. This limiting value will be at the critical 
value when T* = 0, i.e. 

°Zf e T = 2g 2 (1) = 2C (2) = — . (37) 
Nk B T* yw su 3 v ; 

In contrast, for the finite 2D Bose gas, we expect to see a deviation from the above men- 
tioned low-temperature behaviour. In particular, we expect that the reduced heat capacity 
(Cy I 'NksT*) will fall to following the formation of the condensate. Our plots indicate 
distinct differences between the reduced heat capacities for the Bose and Fermi gases, par- 
ticularly at very low T in the vicinity of and below the condensation temperature Tq- Then, 
Cv,Bose is less than Cy^ eTmi and starts to deviate significantly around T c , before rapidly 
going to zero as we approach T = 0. This is in marked contrast to the 3D Bose gas, where 
the conventional heat capacity Cy/Nks peaks at some finite T. In the two-dimensional 
case, however, the heat capacity Cy/Nks is a smooth function of T, in keeping with the 
size-dependent nature of the condensation. Here, instead, we see a "bump" in the plot for 
the reduced heat capacity Cy / NksT* of the 2D Bose gas, indicating the presence of a con- 
densation. Further, as ./V increases, the behaviour tends to that of the 2D Fermi gas as the 
condensation forms at lower temperatures, and consequently, the fall-off of the reduced heat 
capacity occurs at lower temperatures. Indeed, with N = 500000, the two plots of the Fermi 
and Bose gases are virtually indistinguishable. This is significant as conventionally, the 
heat capacity Cy/Nks is used experimentally to determine if a Bose-Einstein condensation 
has occurred. However, as can be seen from above, the heat capacity does not necessarily 



11 



peak for a weak transition which is dependent on the system size, and we have to examine 
alternate expressions to determine the presence of a condensation. 

V. DISCUSSION AND CONCLUSION 

In this paper, we set out to discuss the two-dimensional ideal Bose and Fermi gases in the 
absence of any trapping potentials, and both in the thermodynamic limit and the finite case. 
In particular, we were interested in May's Theorem, and its breakdown for finite systems. 

As can be seen above in Sections 11111 and HV1 the 2D finite Fermi system displays the same 
behaviour as in the thermodynamic limit but the Bose gas deviates significantly for a finite 
system and undergoes a pseudo-Bose-Einstein condensation. It is not a "true" condensation 
as for the 3D Bose gas, where there is a distinct phase transition in the thermodynamic limit 
and a peak in the heat capacity. Rather, here we have a characteristic condensation temper- 
ature Tc that is dependent on the system size and tends to zero in the thermodynamic limit. 
Moreover, though these is no heat capacity peak here, it still deviates significantly from the 
Fermi heat capacity below Tc with the deviations tending to zero as the system size and/or 
the temperature increase. In other words, we have a peak in the reduced heat capacity, as 
shown in Fig. 3. This raises the question of what exactly is defined as a phase transition. 
The two-dimensional Bose gas undergoes a significant change at low temperatures but only 
in the finite case. Thus, it is not what may be termed a "conventional" phase transition 
as for the three-dimensional Bose gas, since it ceases to exist in the thermodynamic limit. 
However, from an experimental perspective, a system is never truly in the thermodynamic 
limit, and so, what we have termed a "pseudo-condensation" is expected. This also, there- 
fore, casts doubt on the use of the heat capacity as a means of identifying phase transitions, 
as discussed above. 

As a final note, we would like to consider some of the outstanding problems. The above 
analysis is a simple one and deals with only the ideal two-dimensional gas. There is still 
considerable debate whether such behaviour is to be found for an interacting two-dimensional 
gas. Some authors^ have claimed that there is a true condensation for the finite interacting 
2D system at sufficiently low temperatures as the coherence length becomes larger than the 
condensate size. Others, however, have claimed that any phase transition would disappear 
in the presence of particle interactions^, as predicted by the Hohenberg theorem. Clearly, 
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the presence of interactions will modify our results greatly. As yet, however, there is no 
indication if there is any correspondence between the two-dimensional Fermi and Bose gases, 
when they are interacting, and further, if a pseudo-condensation is to be found for the finite 
Bose gas. Further, similar analogues exist in other dimensions for ideal gases in a trapping 
potentia l 32 ! 33 . It would be interesting to debate whether there are deeper implications due 
to the Fermi-Bose correspondence, and where such an situation occurs, it is a special unique 
case. This is an area of active work, with new experimental innovations making the resolution 
of the above problems and the elucidation of two-dimensional gases only a matter of time 
now. 
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FIG. 1: Plot of reduced characteristic condensation temperature ksTc/ec against number of 
particles N in the system. [N.B. Spin degeneracy a = 1] 
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Reduced Temperature 

FIG. 2: Plot of condensate fraction Nq/N against reduced temperature ksT jec for (from right to 
left) N = 10 3 , 10 4 , 10 5 , 10 6 , 10 9 , 10 50 . [a = 1] 
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FIG. 3: Plot of difference in the reduced heat capacities Cy/NksT* of the Fermi and Bose gases 
against reduced temperature T* = ksT/ec for (from top to bottom) N = : 500, 5000, 50000, 
500000. [a = 1] 
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